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Abstract 

Attempts at modeling recent expenments of Cheng et al. > indicated that discrepancies 
between theory and experiment can be a result of the form of assumed PDF and/or the turbulence 
model employed. 2 In this tnvestigation, improvements in both the form of the assumed PDF and 
the turbulence model are presented. The results are again used to compare with measurements. 

Initial comparisons are encouraging. 

Introduction 

The successful destgn of scramjet engines requires the understanding of turbulent mixing 
and its effect on combustion at supersonic speeds. The availability of detailed measurements of 
compos, don and temperature and their variances in H 2 -Air supersonic reacting shear layers 
provides the modeler with useful data to develop the codes that will eventually be used in the 
design. The first attempt a, predicting the data of Ref. 1 indicated needed improvements in the 

form of the assumed PDF and the turbulence model. 

The joint-Beta PDF used in Ref. 2 was developed by Girimaji 2 to study mixing. 

Accord, ng ,0 the model, all vanances and covariances tnvolving mass fracnons depend on them 
mean values and an additional quantity which is the sum of their variances. Thus, using this PDF 
requires one additional model equation. The newly developed PDF that is presented below 
results in independent vanances and covariances thus necessttating additional model equations. 
Because of the above, the present effort limits the Joint-Beta PDF to the elemental mass fractions 
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and employs a partial equilibrium assumption to relate the minor species to the elemental mass 
fractions. To accomplish this, an approach similar to that employed by Janicka and Kollmaiuf 
was used in conjunction with Table 1 (Ref. 5) to determine combined variables whose production 
rates are tndependent of the two-body reactions (the fust fourieen reactions in Table 1). The 
minor species, O, H, OH, H 2 0 2 and H0 2 are then expressed in tenns of the elemental mass 
fractions and one of the combined variables using an equilibrium assumption. 

The New Joint-Beta PDF 


The elemental mass fractions obey the relation 


s = x 1 +*2 + x 3 


+ Xt = 1 


( 1 ) 


where Xj is the mass 


fraction of element 'j*. For H 2 -air combustion, the PDF is assumed to have 


the form (reflecting the three elements H 2 , 0 2 , and N 2 ) 


F = i [ X N (1 _ x 3 ) a ^ +Y2) x*' 1 xf 1 8(1 - s)] 

where 8 is the Dirac delta function and c is a constant determined from a normalization 
condition. The parameters a, P, Y„ and y 2 . can be expressed in terms of means, variances, and 
covariances by taking appropriate moments of the joint PDF, i.e., 

a Yi 


( 2 ) 


h = Jxjfdxj 
x? = 

1 


dx 2 dx 3 = 


a + P Yi + Y 2 


- x"j") 2 f dxj dx 2 dx 3 = x^ - (x7) 


x' x ' 2 = I (x 1 -^)(x 2 -x 2 )fdx 1 dx 2 dx 3 = x,x 2 - x t x 2 


( 3 ) 


- a(a + 1) Titti * 1> 

K| “ (a + P)(a + P + 1) (Yi + Y 2 XYi + Y 2 + » 


with 
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g(a+ 1) TO (4) 

XlX2 " (g + fJ)(g + (3+ 1) ’ (Yi+Y 2 TO+Y 2 +1 ) 

— ~2 ~2 ~ 

The parameters g, (3, Yj and y 2 foUow from the expressions for Xj + x 2 , x x , ^ • 


The Combined Variables 


A general procedure has been developed which makes il possible to take a set of 
reactions, such as those indicated in Table 1, and determine the combined variables whose 
production rates are independent of the rates of the two-body reactions. For the set shown tn 

Table 1, the combined variables are 


* 



1 Y h \ Yqh I Yh °2 

2 W H + 2 W 0H 2 W hc>2 


(5) 


Y «2 Y «2 . 
W « 2 = w h 2 



H 


2 W 


H 


1 Yqh _ I Y H°2 _ Y «2°2 

2 W 0H 2 W H02 W H202 


Y H 2 0 Y H 2 0 Io_Ih + I^ 2 +2 
Wh 2 o " W H 2 0 ’ W o W H W H 02 


Y H 2 Q2 

W H 2 02 


( 6 ) 

(7) 


where W is the molecular weight 

The production rates of the combined variables are given by 

= — W 0 2 [Wj 5 + W 16 + W 17 - W 18 ] 
Wh 2 = - 2W H2 [w 15 + w 16 + w 17 - W 18 ] 
W H2 o = - 2W H2 o [W 15 + W 16 + W n - W 18 ] 


( 10 ) 
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The Proposed Model 


The elemental mass fractions. Z Q and Z H and Y^. am dte thtee independent qurntdoes 

solved. The remaining mass fractions, Y 0 , Y H , 


(ID 


for which conservation of mass equations are 
Y OH . Y ho ,. and Y„ 202 , are determined using an equilibrium assumpdon ,n terms of Zq. Z„, 

and y' H 20 . The joint PDF describing temperature and concentration fluctuations ts chosen as 
G(T, Yhj 0 . Zo- Z H- Z N> * f l<D f 2 <Y H 2 0' F(Z » Z H- Z N> 

In the present approach. f,(T) is a Gaussian PDF. f 2 is a single Beta PDF while F is the joint 
Beta PDF indicated in equation (2). 

Results and Discussion 


A schematic of the experimental setup is shown in Ftgure 1 . and Table 2 gives a summ*y 
of tire burner ext. conditions. Measurements were taken « axial locations of 0.85, 10.8. 32.3. 

43 1 , 64.7, and 86. 1 inner diameters downstream of the burner. 

A complex model of the type described above has to be developed step by step. Imrnti 
results of a full finite tare kinetic model calculation with those 
employing a partial equilibrium assumption. These calculations were earned out using a 91 x 91 
grid on a domain of 4 x 4 inchea As a result, only comparisons for the firs, five axiri stations are 
given Figures 2-4 compare measured and calculated mole fractions of HjO and OH and the 
temperature distnbution. As is seen from tire figures, tire partial equilibrium assumption ts an 

acceptable assumption especially a, tire downstream stations. 

The paper will implement the fuU model and will present detailed comparisons wtti, 


calculations compared the 


experiments of Ref. 1 . 
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Tables 

Tab le 1. Chemistry Model 
II Reaction 1 Reaction - 


Number 

—— 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 

16 

17 

18 


1 1 n its nf A 


H 7 + Ch- 
H + Oi - 
OH+H1■ 
O + Hi - 

OH + OH 


OH + OH 
- OH + O 

- HiO 4 - H 
* OH + H 

— HiO 4 - O 


HOi + H — Hi + Oi 
HOi+ H -OH + OH 

HOi + O — Oi + OH 
HOi + OH — HiO + Oi 
HOi 4 - Hi — HiOi 4 - H 
HOi 4 - HOi — HiOi 4 - Oi 
H + HiOi — HiO + OH 
0+ HiOi — HOi + OH 
OH + HiOi — HiO + HOi 
H + OH + M — HiO + M 
H + H + to l — Hi + M 
H + Oi + M — HOi + M 

HiOi + M^OH+OH + M 


. _ ..Uinlu - mnip • S 


170E + 14 

0.00 

24157.0 

. 120E + 18 

-0.91 

8310.5 

.220E + 14 

0.00 

2591.8 

.506E + 05 

2.67 

3165.6 

.630E + 13 

0.00 

548.6 

.130E + 14 

0.00 

0.0 

150E + 15 

0.00 

503.3 

.200E + 14 

0.00 

0.0 

.200E + 14 

0.00 

0.0 

.301E + 12 

0.00 

9411.2 

.200E + 13 

0.00 

0.0 

.100E + 14 

0.00 

1801.7 

.280E + 14 

0.00 

3220.9 

.700E 4- 13 

0.00 

722.2 

221E 4- 23 

-2.00 

0.0 

.730E 4- 18 

-1.00 

0.0 

230E 4- 19 

-1.00 

0.0 

.121E 4- 18 

0.00 

22898.8 
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Figure 3: Mean OH Mole Fraction Comparisons Between Full Finite Rate Chemistry and 
Assuming Partial Equilibrium 




